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ABSTRACT 

We investigate the formation of carbon-enhanced metal-poor (CEMP) stars via the 
scenario of mass transfer from a carbon-rich asymptotic giant branch (AGB) primary 
to a low-mass companion in a binary system. We explore the extent to which material 
accreted from a companion star becomes mixed with that of the recipient, focusing 
on the effects of thermohaline mixing and gravitational settling. We have created a 
new set of asymptotic giant branch models in order to determine what the composi- 
tion of material being accreted in these systems will be. We then model a range of 
CEMP systems by evolving a grid of models of low-mass stars, varying the amount 
of material accreted by the star (to mimic systems with different separations) and 
also the composition of the accreted material (to mimic accretion from primaries of 
different mass). We find that with thermohaline mixing alone, the accreted material 
can become mixed with between 16 and 88% of the pristine stellar material of the 
accretor, depending on the mass accreted and the composition of the material. If we 
include the effects of gravitational settling, we find that thermohaline mixing can be 
inhibited and, in the case that only a small quantity of material is accreted, can be 
suppressed almost completely. 

Key words: stars: evolution, stars: AGB and post-AGB, stars: carbon, binaries: 
general 
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J^ ■ 1 INTRODUCTION tern was an asymptotic giant branch (AGB) star which be- 
came carbon-rich through the action of third dredge-ujjj and 
Carbon-enhanced, metal-poor (CEMP) stars arc denned as r , . , ,, , , , 

. , ,„ ,,-, -m , „ J^"^ „ X,, . — T7~n i „ k . , transferred material on to the low-mass secondary most 

stars with C/FePL>-|-1.0 (IBeers fc Chr isthcb 2005 , with ... . . . ,. . ,. _, . , ,.f v , , 

L ' " ' — ; : — . likely via a stellar wind). 1 he primary became a white dwarl 

Fe/H < —2 m most cases. These objects appear with in- , . , „ , , r . ... ,, , . , 

1 ' ' f -X7*- 11 -1 and has long since taded from view, with the carbon-rich sec- 

creasing frequency at low mctallicity (Lucatcllo ct al. 2006). , , . ,, , . ., , r ,, 

._,_.„_ , : , ondary now being the only visible component ol the system. 

The study of CEMP stars is being used to probe conditions T , , , , . ,. . , . . , . . , 

It has commonly been assumed that the accreted material 
in the early universe. For example, CEMP stars have been . ,, r , ,, , ... ., 

J remains on the surface of the secondary until the star as- 
used to infer the initial mass function in the early Galaxy , , , . , , , , . , . , , , , . r ,, 

. 17 ni^^n ^,1 ■ n , , i- cends the giant branch, at which point the deepening of the 

(e.g. ILucatcllo ct al. 2005a). Chemical abundance studies ,. , , „ , . n , , , , 

; ° : — ■ — ; ~ .' . . convective envelope (referred to as hrst dredge-up because 

have revealed that the majority of the CEMPs are rich in s- . , , , . , . , „,. ,. ... , ,, 

J f- r— 1| 1 material that has experienced CJN-cyclmg in the stellar m- 

process elements like barium (lAoki et al. l 2003). forming the ... .. ,, r s . " ,, . , ... ,, 

' ~ ° tenor is brought to the surface) mixes the matenaf with the 

so-called CEMP-s group. Recent survey work has detected • , • r , 1 TT , 1 r 1 • 1 1 1 1 

. . , .„.. „ , „ T _,, ,_ interior of the star. However, the transferred material should 



a binary companion in around 68% of these CEMP-s stars 



become mixed with the interior of the ac creting star via 



and this is consistent with them all being in binary systems ,, r ,, , ,. . . , rrr — ,.„ ~r\ 

IT 7i n u n n K the process of thermohaline mixing (sec e.g. btancfiftc ct al. 

^Lucatello et al. || 2005b| ). MlChen fc Ha n 2004: BitzaraMet al.ll2003h . This occuTs 

Binary systems provide a natural explanation tor these , , , , , . , , c , , , ,, 

J J ^ r when the mean mofecular weight of the stellar gas increases 

objects, which are of too low a luminosity to have been , ,, c . , , ,. , , , , , 

, , , . „ toward the surface. A gas element displaced downwards and 
able to produce their own carbon. The primary of the sys- 

2 Third dredge-up occurs when the convective envelope of an 

* E-mail: Richard.Stancliffe@sci.monash.edu.au AGB star deepens after a thermal pulse and material that has 

1 [A/B] = log(N\/NB) — log(A r A/A r B)0 experienced nuclear burning is brought to the surface. 
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compressed will be hotter than its surroundings. It will 
therefore lose heat, become denser and continue to sink. 
This leads to mixing on thermal ti mescales until the molecu - 
lar weight difference is eliminated f.Kippenhahn et al.lll980h . 
IStancliffe et al.l (|2007l ) showed that the inclusion of thermo- 
haline mixing could result in the accreted material being 
mixed throughout 90% of the star. 

Recent work has questioned the efficiency of thermo- 
haline mixing in carbon-enhanced meta l-poor stars . Usin, 
a sample of barium-rich CEMP stars, lAoki et al.l <|200; 
showed that the distribution of [C/H] values in turn-off stars 
(i.e. those stars that have reached the end of their main- 
sequence lives, are still of low-luminosity and have yet to 
become giants) was different from that in giants suggesting 
that significant mixing only happened at first dredge-up. 
A simi lar point was made by Dcnissc nkov fc Pinsonneaultl 
(J2008bh using the data of lLucatello et al.l l|2006r i. These au- 
thors showed that the turn-off stars and giants were consis- 
tent with coming from the same distribution if first dredge- 
up resulted in the [C/H] value (and also the [N/H] value) 
being reduced by around 0.4 dex. They find that this re- 
sult is consistent with having an accreted layer of material 
mixed to an average depth of about 0.2 Mq (or alternatively 
having an accreted layer of 0.2 Mq that remains unmixed). 
Neither of these scenarios is consisten t with the extensive 
mixing found bv lStancliffe et al.l (|2007n . 

A possible source for re duced thermohaline mix ing effi- 
ciency has been suggested bv lThompson et al.l IJ2008I ) . These 
authors suggest that the action of gravitational settling will 
alter the composition gradient of the accreting star near its 
surface. Helium will settle from the surface, reducing the 
mean molecular weight at the surface but leading to an in- 
crease in the layers beneath. This produces a small region 
in which the mean molecular weight, /i, decreases outwards 
toward the stellar surface - a situation which is stable to 
thermohaline mixing. This stabilising composition gradient 
(a so-called '/z-barrier') can inhibit the process of thermoha- 
line mixing. 

This paper extends the work of IStancliffe et al,l |2007), 
examining the effect of varying the composition of the ac- 
creted material (mimicking accretion from different masses 
of companion) and the amount of material that is accreted. 
We also investigate the effect that gravitational settling has 
on the extent to which material is mixed. 



2 THE STELLAR EVOLUTION CODE 

Calculations in this work have been carried out using a mod- 
ified version of the STARS ste llar evolution code originally 
developed by lEggletonl ll 19711 ) and updated by many au- 
thors fe.g. lPols et al.lll995T ). Th e version used h e re inc ludes 
the nucle o synth esis routines of IStancliffe et al.l i|2005r ) and 
IStancliffe! (J2005J) , which follow the nucleosynthesis of 40 iso- 
topes from D to 32 S and important iron group elements. 
The c ode uses the opacity routines of Eldridge fc Touti 
2004) , which employ interpolation in the OPAL tables 



Initial No. of Final Core 
mass TPs mass mass 

(Mq) (Mq) (Mq) 



[C/H] [N/H] » 



Iglesias fc Roge rs 1996) and which account for the variation 



1 


13 


0.650 


0.636 


0.661 


-0.643 


0.607 


1.5 


10 


0.892 


0.622 


0.613 


-1.527 


0.621 


2 


11 


1.237 


0.659 


0.648 


-1.611 


0.626 


2.5 


13 


2.422 


0.734 


0.501 


-1.710 


0.619 


3 


16 


1.942 


0.825 


0.264 


-0.036 


0.611 


3.5 


23 


1.871 


0.850 


-0.030 


0.831 


0.628 



Table 1. Details of the final state of the AGB models. The 
columns are: initial mass in solar masses, number of thermal 
pulses the model was evolved through (note that not all the en- 
velope has been removed by the end of the run - see the main 
text for more details), the final mass of the model in solar masses, 
the final H-exhausted core mass of the model in solar masses, the 
final [C/H] value, the final [N/H] value and the mean molecular 
weight, fi, of the ejected material, assuming it is fully ionised. 
The values of [C/H] and [N/H] for the average composition of the 
ejected material will be lower. 



in opacity as the C and O content of the material varies. 

At low temperatures and in the carbon-rich atmo- 
spheres of low-mass asymptotic giant branch (AGB ) stars, 
molecular opacities become important (|Marigdl2002l ). These 



are not included in the lEldridge fc Touu (J200J) opacity rou- 
tines. There are no tables of molecular opacities across a 
range of metallicities currently available (t hough this is be- 
ginning to change, see lCristallo et al.ll2007l. for example), so 
we adopt the procedures of iMarigol (J2002J) to account for 
molecular opacity. Briefly, this involves the computation of 
the abundances of the molecules H2, H2O, OH, CO, CN and 
C2 in the equation of state and adding their contribution on 
to the opacity taken from the regular opacity tables. The 
form of the op acity for the m olecules H2, H2O, OH and CO 
are t aken from Keelev II 19701 ) . while for CN the polynomial 
fit of lScalo fc Ulrichl ( 19751 ) is used. The C2 opacity is as- 
sumed to follow that of the CN opacity. 



3 AGB MODELS 

We have evolved a set of models with masses of 1, 1.5, 2, 
2.5, 3 and 3.5 Mq from the pre-main-sequence without the 
use of convective overshooting. A mixing length of a = 2.0 
has been employed throughout. This value is chosen based 
on calibration to a solar model. All the models were evolved 
using 999 mesh points. The metallicity was set to Z — 10 -4 
and the initial compositi on of the stellar material set to a 
solar-scaled composition IjAnders fc Grevesselll989r ). giving 
the models an iron-to-hydrogen abundance of [Fe/H]= —2.3. 
Thermohaline mixing was included during the evolution and 
the mixing rate was e nhanc ed by a factor of 100 above 
the Kippenhahn et al. (1980) value based on the results of 
ICharbonnel fc Zahnl (2007J), who showed that this enhanced 
rate could reproduce the observed behaviour of certain abun- 
dance trends on the first giant branch. On the asymptotic 
giant branch, the AG B specific mesh spacing function of 
IStancliffe et al.l (|200J) is employed, as is the mixing scheme 
used by these a uthors. P r ior to the AGB, we employ the 
mass-loss rate of Reime rsl (Il975l) wit h r\ = 0.4 while on the 
AGB we use the I Vassiliadis fc Wood] (|l993l ) formula. 

All of the models experience third dredge-up (TDUP). 
The dredge-up efficiency A (which is defined as the ratio 
of the amount of material dredged up to the growth of the 
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core in the preceding interpulse period) is about 0.15 in the 
IMq model, rising to about 0.68 in the 1.5 Mq model and 
is close to unity in the other models. The final outcome of 
the models is given in Table [T] Only the 1 Mq model was 
evolved to the white dwarf cooling track. The remaining 
models all suffered from numerical problems in the stellar 
envelope during the superwind phase, with the exception of 
the 2.5 M@ one which suffered an unrelated numerical prob- 
lem. In each case there was a rapid increase in the stellar 
radius prior to the failure of the model, typically just after 
th e peak of a thermal pu l se. Sim ilar behaviour was reported 
bv lKarakas fc Lattanziol i|2007l ). who suggest that the prob- 
lem may be related to the input physics and particularly to 
the opacities. Further study of this phenomenon is clearly 
warranted. 

With the exception of the 2.5 Mq model, all of the mod- 
els have lost over half their stellar mass and they have all 
entered the superwind phase. We therefore estimate the con- 
tribution of the remaining envelope to the total yield by 
assuming that it is stripped from the star without further 
change to its composition. As mass loss is very rapid in the 
superwind phase we believe that we have missed only a few 
thermal pulses (and associated episodes of third dredge-up) 
at most. We record the final yield, the mass removed and 
the age of the final model for each of the masses considered. 
These details are needed to model the accretion on to the 
secondary as described below. 

The average composition of material ejected from our 
models is displayed in Figure [T] Full details of the total 
mass of each isotope ejected are presented in Table IA1I in 
the appendix. The trends are as one might expect for AGB 
stars. Below 3Mq, the models show an increase in the 12 C 
abundance due to the action of third dredge-up (TDUP). 
For the 3 and 3.5 Mq models there is also an increase in 
the 13 C and 14 N abundances as hot bottom burning sets 
in. Hot bottom burning happens because the convective en- 
velope is deep enough in the star to dip into the top of 
the H-burning shell, allowing CNO cycling to occur during 
the interpulse period. Similarly, there is evidence for the en- 
hanced occurrence of the p-burning Ne-Na and Mg-Al cycles 
in the more massive models. The greatest abundance of 22 Ne 
is seen to occur in models of 1.5-2.5 Mq masses. This iso- 
tope is produced from a-captures on to 14 N in the intershell, 
where it is also destroyed by the neutron-providing reaction 
22 Ne(a, n) 25 Mg. This latter reaction becomes active at tem- 
peratures of around 3 x 10 8 K and tends to favour destruction 
of 22 Ne in the hotter, more massive cores of the higher mass 
models. For a more detailed dis cussion of AGB abundance 
patterns and nucleosynthesis, see lStancliffe fc Jeffervl (|2007l ) 
and references therein. 

We no w compare our mode l s to other s in the litera- 
ture. Both iKarakas fc Lattanziol [|2007h and IHerwid (|2004r ) 
have produced models of the sam e metallicity. The mod- 
els of IKarakas fc Lattanziol ([20071 ') - hereafter KL07 -have 
similar input physics to our own, as they use similar mass- 
loss prescriptions and do not employ the u se of convective 
overshooting. Herwig's models employ the iBlockerl (ll995T ) 
mass loss formula and convective overshooting is employed 
at both the bottom of the convective envelope and the pulse 
driven convection zone. 

KL07 have produced models of 2, 2.5, and 3Mq, along 
with several others that do not have the same mass as ours 
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Figure 2. Comparison of the average composition for isotopes 
in three different evolution codes. The upper panel is for a 
2Mq model; the lower is for a 3Mq model. Triangles indi- 
cated the results of this w ork, squares are from the work of 
Karakas & Lattanzio (2007) - KL07 — and circles are the models 
of IHerwid J2004h designated H04. 



and hence are not considered further. IHerwid (|2004f ) - here- 
after H04 - has produced models of 2 and 3 Mq , along with 
higher mass models which we do not consider here. For the 
two initial masses common to all three works, we have plot- 
ted the average composition of the ejecta for those isotopes 
common to all three codes. The results are displayed in Fig- 
ure [2] We note that there is generally good agreement be- 
tween the 2 Mq models, with few of the isotopes showing 
variations of over one order of magnitude. For many of the 
heavier isotopes, our average compositions tend to lie be- 
tween those of KL07 and H04. The maximum mass dredged 
up after a thermal pulse (TP) is similar to that of H04, but 
we dredge up nearly twice the amount of material as we have 
nearly twice as many pulses with TDUP. Compared to KL07 
we also dredge up a similar amount of material, but their 
model has almost twice as many TPs as ours. The number 
of TPs a model experiences is primarily determined by the 
mass loss that the model experiences. If the superwind phase 
begins earlier, a model will experience fewer thermal pulses. 
The issue of mass loss o n the AGB is a serio us problem for 
stellar models (see Stancliffe fc Jeffervi 2007, for a detailed 
discussion). 

The comparison is less reassuring for the 3Mq models 
(bottom panel of Figure [2}, with a much greater spread in 
abundance predictions. For example, the 14 N abundance dif- 
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Figure 1. The average abundance in the ejecta of the models for selected isotopes from the nucleosynthesis network. The elements on 
the right are the isotopes that are used in the structure code. The isotope 'g' is a sink particle whose abundance relates to the efficiency 
of the s-process. 



fers by nearly 3 orders of magnitude! Again, our models tend 
to lie in between those of KL07 and H04. The patterns of 
the CNO isotopes suggests that the H04 model does not un- 
dergo significant hot bottom burning (HBB) as it has lower 
13 C and 14 N, while those of this work and KL07 do. How- 
ever, the KL07 model has much more 14 N (as well as 23 Na 
from the Ne-Na cycle) than ours due to more mass being 
dredged-up per pulse and there being more thermal pulses 
with TDUP and HBB. The KL07 model typically dredges up 
around 7 x 1CP 3 Mq per pulse (although this value reaches 
1.1 x 1CP 2 for one pulse), compared to around 5 x IO -3 Mq 
for our model. We also note the KL07 model has around 40 
pulses with TDUP, compared to the 13 in our model. 



At the low mass end, we have no direct models to com- 
pare to KL07 as we have produced 1 and 1.5 Mq models, 
while they have made 1.25 and 1.75 Mq models. However, 
we note similar trends. In the lower mass models, the dredge- 
up efficiency drops noticeably and sets of models dredge-up 
similar masses for each episode of TDUP. We note that once 



again, the KL07 models have more episodes of TDUP than 
ours. 



3.1 Cessation of third dredge-up 

The 1 Mq model displays some behaviour not found in the 
other models, namely that TDUP operates only between 
pulses 3 and 8 and ceases until the final pulse, where it oc- 
curs once more. We believe that the reason for the temporary 
cessation of TDUP is a consequence of the changing metallic- 
ity of the envelope. As dredge-up occurs, carbon is brought 
into the envelope making it more metal-rich. It is well es- 
tablished that stellar models of higher metallicity experience 
less efficient third dredge-up which tends to occur only at 
hig her core masses and for more violent thermal pulses (see 
e.g. lKarakas et al.ll2002r ). In the 1 Mq model presented here, 
the envelope metallicity (in terms of the CNO elements) is 
raised up to a sufficiently high level that further thermal 
pulses are too weak to cause the occurrence of TDUP, un- 
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til the last pulse which is significantly more violent than its 
predecessors. 

We have also evolved a 0.9 Mo model in the same way 
as described above. We find that it has an epsiode of third 
dredge-up on its 4 th thermal pulse and then TDUP ceases 
until the 10 th (and final) pulse when TDUP happens once 
more. We note that the pulse strengths for the last pulse 
is logLgr/L Q = 7.84 compared to IogLgf/L© = 7.66 for 
the 4 th pulse. Pulses 8 and 9 have pulse strengths greater 
than that for pulse 4, but no TDUP occurs. This seems to 
support the above hypothesis. 

It should be noted that most authors do not find 
TDUP at the low envelope masses noted her e (see e.g. 
Straniero et al. 20031 ; Karakas fc Lattanziol 120071 ). However, 



Stancliffe fc Jeffervl (J2007J) have found dredge- up at low en- 



velope mass and advanced an argument why this is physi- 
cally reasonable. The efficiency and occurrence of TDUP is 
a contentiou s issue and has been for some time. However, 
IWoodl l|1981f ) has shown that the minimum mass for dredge- 
up to occur is strongly dependent on the mixing length 
parameter, a. He shows that for models with metallicity 
Z = 0.001, the minimum stellar mass for TDUP to occur 
changes from 3.4 Mq when a — 0.67 to just 1.13 Mq when 
a — 1.5. We note that our a is larger than the value of 1.75 
employed by KL07 and this could help to explain the dis- 
crepancy between computations from different codes. While 
this is not proof of the reason why we obtain TDUP at low 
envelope mass, it certainly presents an interesting avenue of 
further study. 



4 CEMP MODELS 

We now turn to the modelling of the secondaries in these 
putative AGB binaries. We assume we are free to choose 
the primary mass, secondary mass and amount of material 
accreted to produce a secondary which has the appropriate 
mass (about 0.8 Mq, see below). The primary mass deter- 
mines the composition of the accreted material and the age 
at which mass is transferred to the secondary. In the case 
of a 1 M© primary, material will be transferred at an age of 
around 5.8 x 10 9 years, while for a 3.5 Mq primary matter 
transfer takes place at an age of around 1.9 x 10 8 years. 
Variations in the amount of mass the secondary accretes 
would presumably be ca used by variations in th e separation 
of the system. Following IStancliffe et alj (|2007l ). we assume 
that matter is accreted from the AGB wind (rather than 
from Roche lobe overflow). In all cases we accrete matter 
at a rate of 10 _6 MQyr _1 which is approximately the rate 
at which the se condary would accret e matter via the Bondi- 
Hoyle process (IBondi fc Hovlelll944r ) from the AGB super- 
wind. 

We evolve a grid of models of stars of 0.7, 0.75, 0.78, 
0.79, 0.795, 0.798 and 0.799 M Q each accreting enough ma- 
terial to make the final mass of the star equal to 0.8 Mq, 
which is approximately the turn-off mass for the halo at the 
current age of the universe. We accrete material from each 
of the AGB models described above on to these stars, using 
the average composition of the ejecta and starting the mass 
transfer at the final age of the AGB model. We evolve two 
separate sets of models: one does not include thermohaline 
mixing, while the other does. 
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Figure 3. Mean molecular weight of the stellar material as a func- 
tion of the fractional mass (the mass above the considered layer) 
for the model accreting 0.1 Mq of material from a 1.5 Mq com- 
panion with thermohaline mixing included. The displayed mean 
molecular weight profiles are: just prior to the onset of accretion 
(solid line), just after accretion has ended (long-dashed line), at 
the end of the main sequence (short-dashed line) and just after 
the end of first dredge- up (dotted line). The total stellar mass 
has increased from 0.7 Mq to 0.8 Mq by the end of the accretion 
phase. 



At the metallicity considered here (i.e. Z — 10 -4 ) a star 
of around 0.7 Mq has only a shallow convective envelope. 
Any material transferred to the star does not become mixed 
into the stellar interior via convection. If a deep convection 
zone did exist (as it does in such stars at higher metallicity) , 
accreted material would rapidly be mixed into the stellar 
interior by convective motions. The impact of thermohaline 
mixing would then be reduced as the mean molecular weight 
of the accreted material would decrease as it becomes diluted 
with the pristine stellar material in the envelope. However, 
at low metallicity no such convective mixing happens and 
thermohaline mixing may efficiently mix the accreted matter 
with the stellar interior. 

To better understand the physical processes involved, it 
is instructive to consider the structural changes that occur 
during the evolution of one of these models. Figure [3] shows 
the mean molecular weight profile at varies points in the 
evolution of the model which accretes 0.1 Mq from a 1.5 Mq 
companion with thermohaline mixing taken into account. 
Initially, the star burns hydrogen in its core, raising the 
mean molecular weight in the central regions while the outer 
regions retail their initial composition (the solid line of Fig- 
ure |3J. The star then accretes material from its companion 
and this material has a higher mean molecular weight than 
that of the original stellar material (long-dashed line). This 
situation is unstable to thermohaline mixing and the mate- 
rial mixes with the lower stellar layers, reducing the mean 
molecular weight of the surface layers and increasing that 
in the interior (short-dashed line). As this happens the star 
continues to burn hydrogen in the core with the mean molec- 
ular weight in this region continuing to rise. Eventually, the 
star evolves to the giant branch and a deep convective en- 
velope develops. This is first dredge-up and it results in the 
mean molecular weight of the surface being homogenised 
over the whole of the convective region. The mean molec- 
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Mi 






M acc 


(M ) 








(M ) 


0.1 


0.05 


0.02 


0.01 


0.005 


0.002 


0.001 


1 


0.37 


0.45 


0.54 


0.57 


0.61 


0.65 


0.67 


1.5 


0.20 


0.28 


0.40 


0.47 


0.51 


0.58 


0.63 


2 


0.18 


0.22 


0.33 


0.41 


0.47 


0.55 


0.61 


2.5 


0.12 


0.22 


0.35 


0.43 


0.47 


0.57 


0.62 


3 


0.15 


0.23 


0.35 


0.42 


0.50 


0.57 


0.62 


3.5 


0.07 


0.17 


0.27 


0.35 


0.43 


0.54 


0.60 



Table 2. Mass coordinate of the greatest extent of thermohalinc 
mixing (in solar masses), as a function of companion mass, Mi, 
and the amount of material accreted, M acc . At first dredge- up, 
the base of the convective envelope reaches a mass co-ordinate of 
around 0.35 Mq at its maximum depth. 



ular weight of the envelope is raised because material that 
has experienced hydrogen burning is brought to the surface. 

Table [2] shows the depth to which thermohaline mixing 
is able to mix the accreted stellar material. The most ex- 
tensive mixing in the above grid mixes accreted matter with 
around 88% of the accreting s tar. This is consistent with the 
work of lStancliffe et al.l (120071 ') , whose model involved accre- 
tion from a highly C-enriched 2 Mq companion. In this work, 
the closest corresponding model does not mix as deeply be- 
cause the companion is not as He- and C-enriched as their 
model was. Only 16% of the star is mixed in the case with the 
shallowest mixing. The depth of mixing increases with the 
amount of material that is accreted. This is to be expected 
as the more accreted material there is, the more pristine 
matter must be mixed with it to reduce its mean molecular 
weight to something comparable to its surroundings. In ad- 
dition, the lower the mass of the primary star from which 
material is accreted, the less efficient the mixing will be. 
This is because less massive stars take longer to evolve, so 
the accreting star has more time to burn its own material, 
raising the mean molecular weight of its innermost regions 
and preventing mixing reaching these depths. 

Figure [4] displays the evolution of [C/H] with luminos- 
ity when accreting 0.001, 0.01 and 0.1 Mo of material from a 
1.5 M© companion. Note that the three models without ther- 
mohaline mixing have different [C/H] values because of the 
existence of a thin convective envelope, which mixes some 
of the material during accretion. Figure [5] displays the evo- 
lution of [C/H] with luminosity when accreting 0.02 Mq of 
material from companions of different mass. These plots are 
representative of the full data set. 

For the models without thermohaline mixing, the only 
episode of mixing occurs between luminosities of around 
logL/Lg = 0.7 — 1.4. This is due to the occurrence of first 
dredge-up as the star ascends the giant branch and a deep 
convective envelope develops. In some of the models there 
is a small amount of mixing into a surface convection zone 
during the accretion phase. In the models with thermohaline 
mixing there is a sharp drop in the [C/H] values at low lu- 
minosities while the star is on the main sequence. If a large 
quantity of material is accreted then it becomes mixed to 
a depth lower than the convective envelope reaches during 
first dredge-up. In this case, we do not see a decrease in 
[C/H] at first dredge-up. If the accreted material does not 
mix as deeply, when first dredge-up occurs the material be- 
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Figure 4. The evolution of the surface [C/H] ratio as a function 
of luminosity when accreting 0.1 Mq (solid line), 0.01 Mq (dot- 
ted line) and 0.001 Mq (dashed line) of material from a 1.5 Mq 
companion. The mass of all the models after accretion has fin- 
ished is 0.8 Mq. The top panel displays models without thermo- 
haline mixing while the lower one is for models with thermoha- 
line mixing. The tracks begin from the p oint at which ac c retion 
finishes. Crosses denoted the CEMPs of Lucatcllo ct al 



2006) 



2007) 



whi le squares denote the Ba-rich CEMPs of lAoki et al 
and lAoki et al.l ll2008h 



comes further diluted and the [C/H] value falls further. This 
can be seen in lower panel of Figure [4] The model which ac- 
cretes 0.001 Mq (the dashed line in the lower panel of this 
figure) only mixes down to a mass co-ordinate of 0.63 Mq 
whereas the convective envelope reaches a mass co-ordinate 
of around 0.35 Mq at its maximum depth. As first dredge- 
up occurs, we see a marked decrease in the [C/H] value as 
expected. 

Both sets of models have problems explaining the ob- 
served patterns of [C/H] as a function of luminosity. Mod- 
els without thermohaline mixing struggle to populate those 
turn-off stars with low [C/H] values and unless a large quan- 
tity of material is accreted, they suffer too much dilution at 
first dredge-up. Models with thermohaline mixing do not 
produce the highest [C/H] values as the accreted matter is 
mixed with the pristine stellar material too quickly to be ob- 
served. In addition, they suggest that we should be able to 
observe turn-off stars with comparatively low [C/H] values. 
No such objects have so far been detect ed and this is no t 
a selection effect (see the discussion in lAoki et alj 120071 ). 
Both model sets favour accretion from 1.5-2.5 Mq compan- 
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Figure 5. The evolution of the surface [C/H] ratio as a func- 
tion of luminosity when accreting 0.02 Mg of material for a 1 Mg 
(solid line), 1.5M (dotted line), 2M (short-dashed line), 3M 
(long-dashed line) and a 3.5 Mq (dot-dashed line) companion. 
The mass of all the models after accretion has finished is 0.8 Mq. 
The top panel displays models without thermohaline mixing while 
the lower one is for models with thermohaline mixing. The tracks 
begin from the point at which accretio n finishes. Crosses denoted 
the CEMPs of iLucatello et all (2006) whi le squares denot e the 
Ba-rich CEMPsof TAoki et al.l l200ik and lAoki et al.l J2008h . 



Figure 6. The evolution of the surface [N/H] ratio as a function 
of luminosity when accreting 0.1 Mq (solid line), 0.01 Mq (dot- 
ted line) and 0.001 Mq (dashed line) of material from a 3.5 Mq 
companion. The mass of all the models after accretion has fin- 
ished is 0.8 Mq. The top panel displays models without thermo- 
haline mixing while the lower one is for models with thermoha- 
line mixing. The tracks begin from the p oint at which ac c retion 
finishes. Crosses denoted the CEMPs of ILucatello et al. (2006) 



whi le squares denote the Ba-rich CEMPs of lAoki et al 
and lAoki et al.l J2008h 



2007) 



ions as accreting different masses from these stars can repro- 
duce the spread in the [C/H] values at luminosities above 
logL/Lg = 1.4. 

The evolution of [N/H] as a function of luminosity when 
accreting a varying amount of mass from a companion of 
3.5 Mq is shown in Figure HI] Figure [7] shows the evolution 
of [N/H] as a function of luminosity when accreting 0.02 Mq 
of material from companions of different mass. The mod- 
els without thermohaline mixing show exactly the same be- 
haviour as was noted for [C/H] , with first dredge-up causing 
a major dilution of the accreted material. In the case that 
the accreted material is nitrogen-poor (as it is for accre- 
tion from companions below 3Mq), there is a slight rise in 
[N/H] at the end of first-dredge-up, following the sharp drop 
caused by the dilution. This occurs because the convective 
envelope reaches down to those regions of the star where 
carbon originally present in the star has been processed to 
nitrogen via the CN cycle. This rise in [N/H] is exactly what 
would be expected in a 'normal' metal-poor star. 

The picture is very different for models including ther- 
mohaline mixing. As with the behaviour of [C/H], there is 



initially a sharp drop in [N/H] on the main sequence as the 
accreted material mixes with the stellar interior. However, 
as the star evolves up the giant branch (i.e. to higher lu- 
minosity), a sharp increase in [N/H] can occur during first 
dredge-up. This occurs because when carbon is mixed deep 
into the stellar interior it can be processed to nitrogen and 
subsequently brought to the surface during first dredge-up 
Iptancliffe et al.ll2007l ). This effect can be seen in the lower 
panel of Figure [7] most notably for the case of accretion 
from a 1.5 Mq companion (the dashed line). The extent of 
this effect depends on how much carbon is present in the 
accreted material and to what depth it is mixed. 

Nitrogen is a real problem for these models. Only the 
3 and 3.5 Mq models produce nitrogen abundances that are 
close to the observed values, with the 3.5 Mq model being 
able to reproduce the spread in [N/H] values. Without ther- 
mohaline mixing, it is difficult to reproduce the low [N/H] 
values in turn-off stars, whereas the thermohaline mixing 
models can populate this region provided that only a small 
amount of material is accreted. It should be noted that 
such massive companions do not reproduce the [C/H] ob- 
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Figure 7. The evolution of the surface [N/H] ratio as a func- 
tion of luminosity when accreting 0.02 MQof material for a 1 Mq 
(solid line), 1.5 M Q (dotted line), 2M (short-dashed line), 3M 
(long-dashed line) and a 3.5 Mq (dot-dashed line) companion. 
The mass of all the models after accretion has finished is 0.8 Mq. 
The top panel displays models without thermohaline mixing while 
the lower one is for models with thermohaline mixing. The tracks 
begin from the point at which accretio n finishes. Crosses denoted 
the CEMPs of iLucatello et alJ 112 006) whi le squares denot e the 
Ba-rich CEMPsof TAoki et al.l l bo07l l and kola et al.l J2008h . 



servations very well. It has been noted that it is hard to 
match b oth the carbon and n itrogen enrichments from AGB 
sources II Johnson et al.ll2007l) . 

The difficulty in matching the models to the observa- 
tions is almost certainly a problem on the theoretical side. 
It is well known that extra-mixing mechanisms must be ac- 
tive during the TP-AGB (and indeed at other phases in the 
lives of stars). Progress at determining the physical cause 
of this unknown extra mixing, dubbed 'cool bottom pro- 
cessing' (CBP) has not been forthcoming. We know that 
some extra-mixing processes must happen in order to get a 
13 C pocket to provide neutrons for the s-process. H owever, 
we no te that our yields are comparable to those of iHerwia 
l|2004l ) who included convective overshooting at the base of 
both the intershell convection zone and the convective enve- 
lope. It is also interesting to note that there is a degree of 
N-richness associated with nearly all CEMP stars. Whatever 
the extra-mixing mechanism is, it may have to produce just 
the right amount of mixing over a range of stellar masses 
which should seriously constrain the physical process. 

One possible alternative explanation for the simulta- 



neous enhancement of both carbon and nitrogen is the oc- 
currence of a proton ingestion episode or flash-driven deep- 
mixing (FDDM). These two labels refer to the situation 
that occurs when protons are ingested into the intershell 
convection zone and burn at very high temperatures. Such 
episodes should prod uce enhancements of both C and N 
IjFuiimoto et al.l |200fj) . In order for the intershell convec- 
tion zone to be able to penetrate into H-rich regions the 
hydrogen burning shell cannot present an effective entropy 
barrier. This only happens at very low metallicity when the 
efficiency of CNO burning is substantially reduced due to 
the absence of CNO nuclei. This does not seem to happen 
at [Fe/H]~ —2 but may account for some of the CEMPs 
with lower metallicity. 

4.1 The effect of gravitational settling 

There are certainly problems with models including thermo- 
haline mixing from an observational perspective. They do 
not seem to reproduce those objects with high [C/H] ratios 
and they imply there should be a population of low [C/H] 
objects that are not observed. However from the theoreti- 
cal perspective there are problems with models that do not 
include thermohaline mixing - one cannot ignore physics! 
If thermohaline mixing does not seem to be effective, there 
mus t be some mechanism t o suppress it. 

iThompson et al,l (|2008l ) have suggested that the effects 
of gravitational settling may inhibit the action of thermo- 
haline mixing. Gravitational settling would lead to helium 
and other heavy isotopes diffusing away from the stellar sur- 
face during the main sequence. This would raise the mean 
molecular weight, fj,, just below the stellar surface, providing 
a positive /i-gradient that would in hibit the action of ther - 
mohaline mixing (see Figure 12 of IThompson et al.l 120081 ). 
However, we note that this /x-barrier is not very extensive in 
their model. At its peak \i is just below 0.6 which is some- 
what less than the [i of the accreted material. For example, 
the mean molecular weight of material ejected by the 1.5 Mq 
model is 0.621. This barrier is also located close to the stel- 
lar surface (in mass). One would therefore expect that the 
accretion of a sufficiently large amount of material will easily 
overcome the barrier. 

To test this, we have implemented the physics of gravi- 
tational settling in the code. The rate of change of the abun- 
dance, Xi, of an isotope i in the absence of nuclear reactions 
and convective mixing is given by 

dXi 1 d . 2 v -. „, 

-W-~^r {rpXiVD) (1) 

where p is the density, r is the radius and vp is the diffusion 
velocity of the species (|Pelletier et al.lll986l ). The diffusion 
velocity is given by 



vd = —Di 



a ln« 



[mi — /i] 



dlnT 



ai2- 



(2) 



dr L n kT dr 

where \i is the mean molecular weight of the material, g is 
the local gravity, a = n(i)/n(H) is the concentration of the 
isotope i relative to hydrogen and rrij is the mass of the 
isotope i (JMichaud fc Va uclair 1991). D12 and 012 are the 
atomic and the r mal d iffusion coefficients and are taken from 
IPaquette et al,l |l98fj). The right-hand side of equation \T\ is 
added to the standard equation for the composition change, 
which is then solved as usual. 
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We re-ran the thermohaline mixing models with grav- 
itational settling included in the calculation. Again, it is 
instructive to consider the structural changes of an individ- 
ual model to understand the processes involved. Figure [8] 
shows the mean molecular weight of the interior of a 0.7 Mq 
star accreting 0.1 Mq of material from a 1.5 Mq compan- 
ion when both thermohaline mixing and gravitational set- 
tling are taken into account. At the zero-age main sequence, 
the composition of the material is uniform as no significant 
H-burning has occurred (the solid line in Figure [8} • Grav- 
itational settling, which acts over long timescales, leads to 
helium (and other heavy elements) settling from the surface 
and hydrogen being enhanced there. This leads to the mean 
molecular weight decreasing in the surface regions (down to 
a fractional mass of around dM/M « 0.01) and a mean 
molecular weight gradient building up (see the long-dashed 
line of the figure) . It is this which may potentially inhibit any 
thermohaline mixing. Note that the mean molecular weight 
remains constant down to around dM/M as 0.01 because of 
the presence of a shallow convection zone. Material of higher 
mean molecular weight is then accreted from the companion 
(short-dashed line). In this example, a large quantity of ma- 
terial is accreted and the mean molecular weight gradient is 
quickly overwhelmed. The accreted material then mixes via 
thermohaline mixing. Once the equilibrium configuration is 
reached, the effects of gravitational settling begin to assert 
themselves again, with helium and heavy elements becom- 
ing depleted in the surface layers by the end of the main 
sequence (dotted line). For this reason, the mean molecular 
weight in the outermost layers is reduced for a second time. 
Finally, the star evolves up the giant branch and first dredge- 
up occurs, homogenising the composition of the envelope 
(dot-dashed line) and producing a similar mean molecular 
weight profile to that of the model without gravitational 
settling. 

Table [3] shows the maximum extent of mixing of ac- 
creted material. We find that gravitational settling is most 
significant when accreting from a low mass companion. This 
is to be expected as settling in low-mass stars acts over long 
(i.e. gigayear) timescales and so it takes time to build up an 
effective mean molecular weight barrier. If accretion hap- 
pens before the barrier is fully established, as it does in the 
case of accretion from all but the two lowest mass compan- 
ions, then thermohaline mixing proceeds in almost the same 
way as it does in the absence of gravitational settling. 

As the amount of material accreted is decreased, the 
effects of the ^-barrier become more pronounced. The ac- 
creted material suffers more dilution with the H-enriched, 
He-depleted surface of the star which reduces its mean 
molecular weight, slowing the rate of thermohaline mixing. 
In the case that the amount of material accreted is very 
small, the barrier can be completely effective and mixing 
does not proceed to great depths. Again, this behaviour 
is expected. The /j-barrier is formed only in the outermost 
parts of the star and if a large quantity of material is ac- 
creted it quickly overwhelms the barrier. When only a small 
amount of material is accreted the barrier is effective at in- 
hibiting thermohaline mixing. If 0.001 Mq of material from 
a 1 Mq companion is accreted, it will only mix with less 
than 1% of the star, in contrast to the 12% seen to occur 
with thermohaline mixing when gravitational settling is not 
considered. 
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Figure 8. Mean molecular weight of the stellar material as a 
function of the fractional mass (the mass above the considered 
layer) for the model accreting 0.1 Mq of material from a 1.5 Mq 
companion with thermohaline mixing and gravitational settling 
included. The displayed mean molecular weight profiles are: at the 
zero-age main sequence (ZAMS, solid line), just prior to the on- 
set of accretion (long-dashed line), just after accretion has ended 
(short-dashed line), at the end of the main sequence (dotted line) 
and just after the end of first dredge-up (dot-dashed line). The 
total stellar mass has increased from 0.7 Mq to 0.8 Mq by the end 
of the accretion phase. The steps in the mean molecular weight 
profile near the surface of the end of main sequence model are 
an artefact caused by the finite precision with which the mass 
co-ordinate is written in the output file. 



Mi 






M acc 


(Mq) 








(Mq) 


0.1 


0.05 


0.02 


0.01 


0.005 


0.002 


0.001 


1 


0.43 


0.55 


0.69 


0.76 


0.78 


0.795 


0.797 


1.5 


0.24 


0.34 


0.46 


0.56 


0.71 


0.795 


0.797 


2 


0.16 


0.25 


0.39 


0.48 


0.56 


0.73 


0.794 


2.5 


0.15 


0.25 


0.39 


0.47 


0.58 


0.74 


0.794 


3 


0.17 


0.26 


0.38 


0.47 


0.58 


0.72 


0.791 


3.5 


0.07 


0.17 


0.30 


0.38 


0.47 


0.64 


0.74 



Table 3. Mass coordinate of the greatest extent of thermohaline 
mixing (in solar masses), as a function of companion mass, Mi, 
and the amount of material accreted, M acc for the models in- 
cluding gravitational settling. At first dredge-up, the base of the 
convective envelope reaches a mass co-ordinate of around 0.35 Mq 
at its maximum depth. 



This effect can also be seen in Figure [9] which displays 
the evolution of [C/H] as function of luminosity when ma- 
terial is accreted from a 1.5 Mq companion. The decline in 
[C/H], followed by the sudden rise around logL/Lg = 0.7 is 
caused by gravitational settling on the main sequence both 
reducing the surface carbon abundance and enhancing that 
of hydrogen. The sudden rise occurs as first dredge-up sets 
in, with the deepening convective envelope mixing the stel- 
lar interior. When accreting 0.1 Mq the material is quickly 
mixed in to the same degree in both models. If less mate- 
rial is accreted, we begin to notice differences in the degree 
of mixing. In the case of accreting 0.001 Mq the /j-barrier 
keeps the [C/H] ratio over 1 dex higher than in the case with 
thermohaline mixing alone. In fact, the /x-barrier keeps the 
[C/H] value higher on the main sequence (logL/Lg < 0.7) 
than it is in the case when 0.01 Mq is accreted. This is in 



10 R. J. Stancliffe & E. Glebbeek 



o - 






-1 - 



V x D □ h d B x a □ □ 


: \ Q \ x- n d x <f n n n x >$ x nx x x 

'A < ' 
•« ? *> 

''A ' V 




- f U « J * D - 


- '■, a >0 n x □ □ □ x ^ x □>; x y 

- i d **b □ * * X ° * ^ * ^ X x - 

- '. x x xx X D * * * 




- 1 xx nx D x - 
1 x 

1 * 

- 1 X 


" , ■ 1 1 " 



1 2 

logj L/L s 



Figure 9. The evolution of the surface [C/H] ratio as a function 
of luminosity when accreting O.lMg (solid line), 0.01 Mq (dot- 
ted line) and 0.001 Mq (dashed line) of material from a 1.5 Mq 
companion. The mass of all the models after accretion has fin- 
ished is 0.8 Mq. The top panel displays models with thermoha- 
line mixing and gravitational settling while the lower one is for 
models with thermohaline mixing only. The tracks begin from the 
point at which a ccreti on finishes. Crosses denoted the CEMPs of 
iLucatello et al.l 1 2006) whi le squares denot e the Ba-rich CEMPs 
of lAoki et alj J2007I) and lAoki et all J2008h . 



contrast to the situation with thermohaline mixing alone, 
where the model which accretes 0.001 Mq (the dashed line 
in Figure^ has a [C/H] value about 0.6 dex lower than that 
of the model which accretes 0.01 Mq throughout its whole 
evolution. 

With the inclusion of gravitational settling, we no longer 
obtain models with low [C/H] values at low luminosities 
(provided the initial [C/H] value is high enough). How- 
ever, gravitational settling does present an added prob- 
lem. As carbon settles from the stellar surface during the 
star's main-sequence lifetime, the [C/H] ratio drops. This 
makes it difficult for us to reproduce the more luminous 
turn-off objects with [C/H] values in the range -1-0. This 
is a serious problem for any model that includes gravita- 
tional settling. If more material is accreted, then a higher 
[C/H] value results throughout the main sequence. For ex- 
ample, if 0.2 Mq of material were to be accreted from 
a 1.5 Mq companion, the resulting object would have a 
[C/H] value with a minimum of just -0.5. We note that we 
have not included the effects of radiative levitation. How- 
ever, based on the simulations of I Richard et al.l (2002) and 



iRichard. Michaud fc Richerl |2002), it seems that radiative 
levitation should be ineffective for carbon. 

There is also still the problem of the high [C/H] ob- 
jects. The highest [C/H] values seem to correspond to those 
models which have accreted the most mass from their com- 
panions. In these cases, the /x-barrier is ineffective and ther- 
mohaline mixing proceeds unhindered. The highest [C/H] 
value obtained in our models is about [C/H] ~ —0.2. If 
more material is accreted, we can reach a higher [C/H]. For 
example, accreting 0.2 Mq of material from a 1.5 Mq com- 
panion produces [C/H] ~ 0. This is in line with the u pper 
[C/H] valu e s of t he Ba-rich stars of lAoki et al.1 (|2007l ) and 
Aoki et al.l (120081) . Ho wever, the sample of C-rich stars of 
Lucatello et all (|2006h has three objects with [C/H] > 0.6. 



It is difficult to explain these objects with models involv- 
ing thermohaline mixing as our AGB models only just pro- 
duce [C/H] values of around this. Any amount of mixing 
will easily reduce this value. However, we note that these 
three high [C/H] CEMPs are all somewhat more metal-rich 
(with [Fe/H] greater than around -1.8) than the models of 
this work so these models may not be directly comparable. 



5 DISCUSSION 



As suggested by iThompson et al.l l|2008l ) a means of sup- 
pressing the action of thermohaline mixing is necessary to 
prevent the destruction of lithium in some CEMPs. Lithium 
is a fragile element that is easily destroyed at temperatures 
of around 2 x 10 6 K which are found at around 0.01 Mq be- 
low the stellar surface of the low-mass secondary while on 
the main sequence. As we have seen, gravitational settling 
will only allow the accreted matter to remain at the sur- 
face if the mass of accreted materi al is small. In the case o f 
CS22964-161, the CEMP binary of IThompson et all [|200SD . 
the measured [C/H] value of -1.2 fits well with accretion from 
a 1-1.5 Mq companion when about 0.001 Mq is transferred 
to the recipient star. We would expect lithium to survive 
in this model because mixing does not reach down to re- 
gions where the temperature is high enough for lithium to 
burn. However, we have not tried to model this. We leave the 
treatment of the behaviour of the light elements to future 
work; 

iDenissenkov fc Pinsonneaultl l|2008bl ) suggested that 
the shift in the [C/H] and [N/H] values at first dredge-up 
requires that the average CEMP has a depth of 0.2 Mq of ac- 
creted material. They point out that this can either be from 
accreting 0.2 Mq of material that remains unmixed, or from 
accreted material mixing to a depth of 0.2 Mq. This would 
imply that, on average, no more than about 0.01 Mq could 
be accreted from a 1 Mq companion (this value rises slightly 
to around 0.05 Mq when gravitational settling is taken into 
account), with this value decreasing rapidly as the compan- 
ion mass increases. This assumes that thermohaline mixing 
is as efficient as is presented in this work. If thermohaline 
mixing is less efficient, then a greater amount of material 
coul d be accret ed. 

lAoki et al.l (|200Sf ) have suggested that the [C/H] val- 
ues measured for turn-off stars represent the composition of 
the donor AGB star. This cannot be reconciled with models 
including thermohaline mixing, as mixing leads to a large 
decrease in the [C/H] value. For the star to retain the com- 
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Figure 10. The distri bution of [C/H] for tu rn-off stars and gi- 
ant sta rs for the da ta of Lucatell o et al.l 11200 6) and the combined 
data of lAoki et all ||2007|) and I Aoki et alf ifeoOS). Stars are iden- 
tified as turn-off objects if logL/LQ < 0.7 (i.e. they have not yet 
undergone first dredge-up) and are displayed in the upper panels. 
Stars are defined as post-first dredge-up giants if log L/Lq > 1.5 
and are displayed in the lower panels. The total number of stars, 
Ntot, in each plot is also displayed. 



position of the AGB donor a substantial amount of material 
would have to have been accreted from the primary, such 
that the accreted material dominates the original stellar ma- 
terial. However, the accretion and subsequent thermohaline 
mixing of a large quantity of material does not agree with 
the observed differences in distribution of [C/H] from turn- 
off to giant stars (see Figure 1101 right-hand column) . Effi- 
cient thermoahline mixing of a large quantity of accreted 
material leads to very little change in [C/H] at first dredge- 
up. This cannot be the case in the majority of CEMP stars. 
In addition, if [C/H] ~ does represent the composition 
of the donor then there is a major problem with the AGB 
models which predict substantially larger values. This is not 
improbable given how little is known about the mass-loss 
rate at low metallicity and the efficiency of third dredge-up. 
The existence of high [C/H] CEMPs suggests that ther- 
mohaline mixing is not efficient in all cases. The difference in 
the distri bution of [C/H] fo r turn-off and giant stars in the 
sample of lAoki et al.l i|2008l ) suggests that some, but not all, 
turn-off stars with [C/H] m do not become mixed after first 
dredge- up while some do (see Figure [lO)l . After first dredge- 
up there is still a large proportion of stars with [C/H] ~ 
which suggests that little dilution of material has taken place 
in these objects during first dredge-up. These objects must 
have experienced efficient thermohaline mixing, while the 
others remained unmixed to some degree until first dredge- 
up. T here is less evid e nce fo r this phenomenon in the data 
set of lLucatello et al.l (|2006l ). which displays a distinct shift 
in the peak [C/H] value between turn-off and giant stars (see 
the left-hand column of Figure HOfl . In order to predict the 
distributions of [C/H] for turn-off stars and giants based on 



the models presented herein, one would need to do a pop- 
ulation synthesis calculation for each model set. Attempts 
at modelling the CEMP population are currently under way 
(Izzard et al., in prep). 

If thermohaline mixing cannot be inhibited by gravi- 
tational settling in the case that a large quantity of ma- 
terial is accreted, we must find another way to reduce its 
efficiency. A possi ble mechanism to do this ha s been put 
forward by iDenissenkov fc Pinsonneaul t (2008a). These au- 
thors point out that rotationally-induced horizontal turbu- 
lent diffusion can suppress thermohaline convection and ex- 
amine the specific case of mixing on the giant br anch (see e.g. 
lEggleton et~ai]|2006| ; [Charbonnel fc Zannll2007h . The molec- 
ular weight inversion found there is at least two orders of 
magnitude smaller than in the case of accretion in CEMPs, 
so one presumes the turbulent diffusion would need to be 
much stronger to suppress thermohaline mixing in this situ- 
ation. The interaction between thermohaline mixing and ro- 
tation is an intriguing possibility and merits further study. 
We are unable to investigate this at present as our code does 
not include rotational physics. 



6 CONCLUSIONS 

We have modelled the accretion of AGB material on to low- 
mass stars in order to determine what chemical signatures 
may be observed in CEMP stars. We have examined three 
specific cases: canonical evolution including only convective 
mixing, the inclusion of thermohaline mixing and the inclu- 
sion of both thermohaline mixing and gravitational settling. 
We find that thermohaline mixing can lead to accreted ma- 
terial being mixed with between 16 and 88% of the accreting 
star, depending on the mass and composition of the accreted 
material. When gravitational settling is included, thermoha- 
line mixing is severely inhibited when only a small amount 
of material (around a few 1O _3 M0) is accreted because of 
the presence of a /z-barrier formed by the settling of helium. 
This barrier is less and less effective as the amount of ac- 
creted material is increased. It is also less effective for more 
massive companions as the barrier has had less time to be 
established before accretion on to the secondary occurs. 

Models without thermohaline mixing produce turn-off 
objects with [C/H] values that are too high to match obser- 
vations. Unless a substantial quantity of material has been 
accreted, these models also suffer from too much dilution at 
first dredge-up. Models with thermohaline mixing cannot re- 
produce the highest [C/H] values observed and predict the 
existence of low [C/H] objects which are not observed. In 
order to reproduce the highest [C/H] values, which seem to 
require the accretion of a large quantity of material from a 
companion, some alternative mechanism to suppress ther- 
mohaline convection is required. The inclusion of gravita- 
tional settling can solve the problem of the low [C/H] objects 
at low luminosity as the ^-barrier prevents mixing when the 
amount of accreted material is small. However, the inclu- 
sion of this physics presents another serious problem namely 
that carbon will settle from the surface during the main se- 
quence, making it extremely difficult to form turn-off stars 
with [C/H] values of -1-0. 

There also appears to be a problem with the abundance 
predictions of the AGB models. They do not predict sub- 
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stantial enhancements of both carbon and nitrogen at the 
same time (except in a very narrow mass range). Further 
work needs to be done to find a mechanism capable of pro- 
ducing the observed abundance patterns. 
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APPENDIX A: STELLAR YIELDS 

Here we present the final gross yields (i.e. the total mass of 
each isotope ejected) from each of the AGB models for all 
of the isotopes in the nucleosynthesis network. 
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Table Al. Final gross yields in stellar masses for all the isotopes in the nucleosynthesis network. 



